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LiCr02 is a 2D triangular antiferromagnet, isostructural with the common battery material 
LiCo02 and a well-known Jahn- Teller antiferromagnet NaNi02. As opposed to the latter, LiCr02 
exibits antiferromagnetic exchange in Cr planes, which has been ascribed to direct Cr-Cr d — d 
overlap. Using LDA and LDA+U first principles calculations I confirm this conjecture and show 
that (a) direct d ~ d overlap is indeed enhanced compared to isostructural Ni and Cr compounds, 

(b) p — d charge transfer gap is also enhanced, thus suppressing the ferromagnetic superexchange, 

(c) the calculated magnetic Elamiltonian maps well onto the nearest neighbors Heisenberg exchange 
model and (d) interplanar inteaction is antiferromagnetic. 



The series of compounds with a common formula 
AAI02, where A is an alkaline or noble metal, usually 
Li or Na, and M is a 3c? metal, formed by triangular 
AIO2 layers stacked hexagonally {e.g., LiCo02) or rhom- 
bohedrally {e.g., LiNi02) with full or partial intercala- 
tion by A, has been attracting considerable recent in- 
terest, largely driven by the immense importance of the 
LiCo02 compound in electrochemical industry and by 
the unconventional superconductivity discovered in the 
hydrated Nai/3Co02. The nickelates, LiNi02, NaNi02, 
AgNi02, Ag2Ni02 have been also subject of numerous 
studies, mostly because of their magnetic properties cou- 
pled with interesting structural transformations. How- 
ever, chromates, such as LiCr02, NaCr02, and KCrOg, 
despite their potential use for rechargeable batteries [1| 
and as catalists 2] have been studied experimentally only 
sporadically^, A, ^, j), _7, .8.], and no first principle calcu- 
lations, to the best of my knowledge, have been reported 
so far. 

In this paper I report all-electron full-potential elec- 
tronic structure calculations for LiCr02. In agreement 
with the experimental findings, the magnetic interac- 
tion in-plane is found to be strongly antiferromagnetic. 
Interplane magnetic interaction is very weak and also 
antiferromagnetic. The total energy calculations for 
three different collinear magnetic configurations map per- 
fectly well onto the standard nearest-neighbor Heisenberg 
model. As conjectured in the first experimental papers[3| 
the main reason for switching the in-plane magnetic in- 
teractions from ferromagnetic in LiNi02 to antiferromag- 
netic in LiCr02 is mainly the enhanced direct overlap 
between the metal d— orbitals in chromates, while, addi- 
tionally, the increased charge-transfer p — d gap reduces 
the ferromagnetic superexchange in chromates as wellQ. 
Finally I will discuss the role of Coulomb correlations as 
revealed by LDA+U calculations. 

LiCrOo crystallizes in a rhombohedral R3m 
structurejJl with the lattice parameters a — 2.898 
A, c = 14.423 A, and with the O height zq = 0.261. 
Its magnetic structure is close to the ideal 120° struc- 
ture characteristic of the nearest neighbor Heisenberg 
model on a triangular plane. The computational re- 



sults reported below were obtained using the standard 
full-potential linearized augmented plane wave code 
WIEN2k[l3|. Exchange and correlation were taken into 
acount in a gradient approximation of Perdew et al [ll| . 
Convergence with respect to both the cutoff parameter 
RKraax and the number of the inequivalent k-points (up 
to 400-1- ) was checked. 

I start first with the hypothetical ferromagnetic struc- 
ture. Cr d-states are split by the crystal field into 3 t2g 
and 2 eg states, separated by roughly 3 eV in the spin- 
majority and 2 eV in the spin- minority channel (Fig. [T]). 
The difference is due to the fact that the crystal field due 
to 0-Cr hybridization is stronger in the spin-majority 
channel where the energy separation between the oxy- 
gen p-states and the Cr d— states is smaller. Cr^+ has 
exactly three electrons, and due to the strong exchange 
splitting it turnes out to be well insulating already in the 
LDA; thus it can be treated as band insulator. (I will 
show later that the effect of the Mott-Hubbard correla- 
tions is relatively small), with the band gap between the 
occupied spin-up and unoccupied spin-down t2g bands. 
The minimum spin-flip gap is slightly larger than 1 eV. 
The average spin-flip energy between the two bands is 
close to 2.9 eV. The oxygen states are well separated 
from the d— states, as opposed to the isostructural oxides 
of the late 3d metals (Co, Ni), where the correspond- 
ing bands essentially overlap. The calculated magnetic 
moment is, obviously 3 /is, not too far from the ex- 
perimentally measured noncoUinear magnetic moment of 
2.68±0.13 /is ^5]. Structural optimization with experi- 
mentally constrainted cell dimensions a and c leads to 
the O position zq — 0.258, corresponding to an approxi- 
mately 4% larger distance between Cr and O planes com- 
pared to the experiment [4] , and a Cr-O-Cr angle of 93° 
compared to the experimental 94.6°. I will discuss later 
possible relation of this deviation to magnetic properties. 

LiCr02, being a good isulator already in LDA, should 
exhibit an antiferromagnetic superexchange interaction 
between the Cr02 layers, proportional to t\ / j. , where 
t_L is the effective interlayer hopping, a small number 
since it proceeds via a long path of Cr-O-O-Cr^, and 
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FIG. 1: Band structure of LiCr02 in the hypothetical fer- 
romagnetic state. Spin-majority bands are solid and blue, 
spin-minority ones dotted and red. (color online). 

As./, is the energy cost for flipping a d electron spin (in 
LSDA, it is Stoner /, in the Hubbard model U). To 
get an estimate (from above [l3j|!) of this interaction I 
compare the total energy of the ferromagnetic state and 
the A-type antiferromagnetic state (ferromagnetism in 
plane, antiferromagnetic stacking). The energy difference 
appears to be 3±1 meV per two formula units, a very 
small number indeed. This strongly suggests that LiCr02 
is a very good model material for studying frustrated 
magnetism in two dimensions. NaCr02 maybe even a 
better example of the same physics, as discussed in Ref. 
Q from an experimental point of view. 

Let me now discuss magnetic interaction in plane. The 
above-described exchange splitting that renders LiCr02 
a band insulator in LDA calculations, being essentially a 
local effect, should be operative in an antiferro-, as well 
as in any ferrimagnetic arrangement. In other words, 
no new magnetic physics is introduced by Mott-Hubbard 
effects. Indeed, in the antiferro- and ferrimagnetic cal- 
culations, described in more details below, I invariably 
found insulating ground states (a corollary of this find- 
ing is that the magnetically disordered state above T/v 
will be also insulating even in LDA). 

The isostructural compounds with higher transition 
metals (Co, Ni), when order magnetically, assume the so- 
called A-type antiferromagnetism, that is, ferromagnetic 
planes, stacked antiferromagnetically. The mechanisms 
for the in-plane ferromagnetism are well understood: this 
is the classical 90° superexchange, plus, in metallic com- 



pounds, like Nao.7Co02, Stoner ferromagnetism. LiCr02 
in not metallic, however, the 0-Cr-O bond angle is fairly 
close to 90° and one expects the corresponding superex- 
change to be ferromagnetic. 

A popular explanation [Sj] is that this superexchange is 
indeed ferromagnetic but is surpassed by direct-overlap 
exchange between Cr d-orbitals. It was also pointed 
out@ that the O p band in the higher metal (Co, Ni) is 
located higher with respect to the metal d-band, which 
enhances the p — d hybridization and the superexchange, 
compared to that in LiCr02. While logical per se, how- 
ever, neither of these propositions has been subjected 
to to a quantitative test. Besides, although among the 
physicist dealing with transition metal oxides it is gen- 
erally known that direct d — d exchange is antiferromag- 
netic, this is not a trivial or obvious statement. Suffice 
it to remind the reader that essentially all textbooks in 
solid state physics describe direct exchange in terms of 
the original Heisenberg theory that in fact leads to the 
ferromagnetic interactions. It is instructive to revisit the 
issue of the net sign of a direct exchange interaction from 
the local spin-density functional (LSDA) point of view, 
which is the foundation of all quantitative investigations 
in this direction, and I refer the reader to the Appendix 
where such an analysis is presented. 

To address the nature of the magnetic interactions in 
LiCr02, let us compare it with the isostructural LiNi02. 
LSDA calculations predict in the latter a ferromagnetic 
exchange in plane, of the order of 5 meV[15l]. To com- 
pare with LiCr02, I have performed LDA calculations in 
a 2x2 X 1 supercell, assuming three different magnetic 
patterns inside the plane: ferromagnetic, antiferromag- 
netic (alternating ferromagnetic chains running along 
100, see Fig. and ferrimagnetic, relatining the hexag- 
onal symmetry by flipping one spin out of four (Fig. 
The results can be perfect well mapped onto the near- 
est neighbor Ising model, yielding antiferromagnetic ex- 
change > 20 meV (22 or 23 meV, depending on which 
two lines in Table 1 are used). This is, of course, in ex- 
cellent agreement with the experiment, but how does it 
answer to the theoretical conjectures described above? 




FIG. 2: Two different magnetic pattern ("antiferro": left, 
"ferrimag: right" ) used in calculating the exchange constants. 
Filled (open) circles indicate up (down) moments within in 
the supercell, solid (dashed) line ferro- (antiferro-)magnetic 
bonds. 

To this end, it is instructive to switch to a less accurate, 
but more flexible method. Linear MT-orbitals (LMTO), 
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TABLE I: LSDA energies with respect to the energy of the 
FM state (in meV), magnetic moments insides the Cr MT 
spheres (in /is), and excitation gaps (in eV) for several mag- 
netic states in a quadruple 2x2 LiCr02 supercell. 
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which provides a capabihty of eliminating particular ra- 
dial functions from the basis, and thus assessing effects 
of particular orbital overlaps on the bands. I have cal- 
culated the band structure of LiCr02 and LiNi02 re- 
moving all oxygen orbitals from the basis, so that the 
resulting d-band width is essentially coming from direct 
d — d overlaps. The resulting band widths appear to 
be for LiCr02 and LiNi02, respectively, 1.4 eV and 0.5 
eV, indicating that the direct d — d overlap integrals for 
Cr are three times large than for Ni, in accord to earlier 
conjectures^!- Comparing these numbers with the band- 
widths in full calculations shows that in LiCr02 most of 
the total bandwidth is due to the direct overlap, while in 
LiNi02 it mainly comes from the indirect hopping via 
oxygen. Note that in some crude approximation this 
means that the direct exchange due to d-d overlap is 
an order of magnitude stronger in LiCr02 compared to 
LiNi02. 

Let us now look at the 0-p ~ metal d energy separation. 
Again, in LMTO there is a gauge that can be readily 
used: separation between the corresponding band-center 
quasiatomic parameters. In LiCr02 this separation ap- 
pears to be 5.5 eV, while in LiNi02 it was 1.6 eV, 3.5 
times smaller. This leads to a suppression (albeit not by 
the same factor) of the ferromagnetic superexchajige in 
the former compound, as proposed by Khomskii 9]. It is 
somewhat hard to estimate the actual reduction of this 
interaction, because the more diffuse character of the Cr 
d orbitals leads to some enhancement of the p — d hopping 
compared to the Ni compounds. Since the p—d overlap is 
already strong in the nickelate, this enhancement is not 
nearly as dramatic as that of the direct d — d hopping, 
but it should be present. 

The summary of this part is that while in the late 3d 
metal layered oxides, such as (Li,Na)Ni02, (Li,Na)Co02, 
the relatively strong ferromagnetic interactions (90° su- 
perexchange, but also Stoner ferromagnetism, discussed 
below) is largely, but not entirely compensated by the 
AFM superexchange (due to deviation from 90°) and by 
the direct d ~ d AFM exchange. In the early 3d metal 
oxides, (Li,Na)Cr02, the latter is greatly (up to an order 
of magnitide) enhanced, while the former is suppressed. 

I will now address the question of the relative impor- 
tance of several factors favoring AFM in-plane interac- 
tions in the early 3d oxides. Two have already been 



mentioned, stronger deviation of the metal-oxugen-metal 
angle from 90° {e.g., in LiCr02 this angle is 94.6°, while 
in NaNi02 is close to 92°), and stronger direct d-orbitlas 
overlap. There is a third important factor: the Stoner 
magnetism. Indeed, in such a system as NaNi02, where 
Ni is a Jahn- Teller ions with an orbital degeneracy, the 
Fermi level falls inside the eg band, which is wider in the 
FM case, and, correspondingly, has lower kinetic energy. 
In other words, itinerant eg electrons have more free- 
dom to move in the crystals on a background of parallel 
magnetic moments (one can also call this a d-d double 
exchange; the difference is purely terminological). Note 
that it does not matter if a small Jahn- Teller gap opens 
up in the eg band; as long as this gap is smaller than or 
comparable to the band width, kinetic energy will still 
favor the FM arrangement. LiCr02 does not display any 
orbital degeneracy, thus lacking this contribution to the 
magnetic interactions. 

It is tempting to try to get an idea of the relative im- 
portance of the two FM interactions. I attemted to ad- 
dress this issue by doing calculations for a hypothetical 
LiCr02 with oxygen octahedra unsqueezed so as to have 
Cr-O-Cr angle exactly 90°, while keeping the intraplane 
Cr-Cr distance constant (otherwise the calculations pro- 
ceeded exactly as for the experimental structure, as de- 
scribed above in Table 1 and Fig. [2|) . One may think that 
this procedure would substantially enhance the 90° FM 
superexchange, hopefully keeping other magnetic inter- 
actions unchanged. The result seems, on the first glance, 
unusual: the exchange interaction becomes more antifer- 
romagnetic (by about 17 meV). A closer look, however, 
reveals that this can be traced down to the fact that 
the equilibrium 0-Cr bond length in the AFM struc- 
ture appears to be longer than in the FM structure (see 
FiglSI, providing less hybridization and slightly smaller 
exchange splitting; note that the AFM energy gain due 
to the Cr-Cr direct exchange inversly depends on the ex- 
change splitting. The total energy curves in Fig. [3] are 
practically rigidly shifted with respect to each other: the 
curvature is the same within the computational accuracy, 
corresponding to the frequency of the Aig phonon of 567 
cm^^. 

Finally, although one can expect that the role of the 
on-site Mott-Hubbard correlations should be small, given 
large widths of the Cr d bands, it is of interest to esti- 
mate it u sing a standard implementation of the LDA-I-U 
technique [lO|[. The latter is most often used in one of the 
two flavors |16|: one is designed to reporoduce the fully 
localized limit (FLL), the other simulates fluctuations 
around the mean- field solution (AMF). The former is be- 
lieved to be more appropriate for large U system, and the 
latter for small J7's (admittedly, in this regime the very 
concept of LDA-I-U becomes questionable). I have esti- 
mated the parameters using the LMTO internal quasi- 
atomic loopfr?! to be [/ = 2.3 eV, J = 0.96 eV. This is a 
moderate C/, probably more on the AMF side. Because of 
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O position in iattice coordinates 
0.258 0.256 0.254 



U=2.3 eV, J=0.96 eV 





FIG. 3: Total energy as a function of oxygen position. The 
upper curve is for the ferromagnetic, the lower for the ferri- 
magnetic case (see Table 1). The range of positions spans the 
experimental structure from Ref.Q on the left to the structure 
with the Cr-O-Cr angle equal to 90° on the right. Calculated 
equilibrium positions are indicated by arrows. 



that, I have performed LDA+U calculations in both lim- 
its, in order to compare them with each other and with 
LDA. The results (assuming the FM structure) are shown 
in Fig. 21 As expected, LDA+U does not introduce any 
new physics: the minimal gap is still the spin-flip gap be- 
tween t2g and eg, enhanced by 0.5 eV in AMF, due to a 
shift of the occupied band by w (C/ — J)/2, and by 1 eV in 
FLL, due to an additional shift of the unoccupied states 
by the same amount. Note that these shifts are smaller 
than those expected in a typical band insulator due to the 
well-known density-derivative discontinuity 18[ , related 
to an unscreened long-range exchange interaction[l^. 
For example, for ZnS, which has a gap of 3.8 eV, LDA 
gives a gap of 1.7-1.8 eV. While I was not able to locate 



data for the experimental gap in LiCr02, the o ptic al gap 
in NaCr02 was reported to be about 3.5 eV[20|, thus 
making comparison with ZnS rather meaningful. Note 
that the numbers in Fig. [4] are for the FM ordering; in- 
troduction of antiferromagnetism of course increases the 
gap (cf. Table 1). 

To conclude, one can describe LiCr02 as a weakly cor- 
related band insulator with the minimal gap beng the 
spin-flip t2g — Cg gap. Magnetically it is a frustrated 
nearest- neighbor 2D triangular antiferromagnet, with a 
sizeable exchange constant of the order of 20 meV. 

I would like to acknowledge many useful discussions 
with Michelle Johannes and Daniel Khomskii. 



APPENDIX: DIRECT EXCHANGE IN LDA 

Consider two overlaping orbitals (jfi = (j){r), <j)2 = 
(/)(|r — R|) centered at the points separated by R. The 
corresponding atomic (on-site) energy we shall call E, 



FIG. 4: Band structure of ferromagnetic LiCr02 in LDA and 
in two different LDA+U flavors, (color online). 



and we assume that there is intraatomic exchange split- 
ting such that the spin-up state has the energy E, but the 
spin-down state has the energy E + I. The intraatomic 
exchange parameter, /, in LSDA is known as "Stoner 
parameter". For strongly correlated systems Hubbard 
U >> I plays the leading role instead. One may ex- 
pect that in an intermediate regime the cost of flipping a 
spin is intermediate between the LSDA I and Hubbard U. 
The so-called Andersen's force theoremfl^ states, among 
other things, that in LSDA the total energy difference be- 
tween the ferro- and antiferromagnetic arrangement can 
be computed as the difference in one-electron energy plus 
the difference in magnetic ("Stoner") energy. The one- 
electron part favors antiferromagnetism. Indeed, the oc- 
cupied levels in the ferromagnetic case are not shifted 
with respect to the on-sites energy for the spin-up state, 
but in the antiferromagnetic case, in the second order 
perturbation theory, they shift down by t^// each, where 
t is the hopping integral, t = (^(pi \ — + Vef f {r)\(j)2) (in 
Ry units). 

The magnetic part favors ferromagnetism, since the 
corresponding energy in LSDA is Est = — J Im'^{r)/A, 
where m^(r) is the total spin density. For the ferromag- 
netic alignment m{r) = (j>l+ (p^, Est = - J I{(j>l+(j>l)'^ /4:. 
For the antiferromagnetic one m{r) = (t)\— (j>2, Est 

-II 

by I \vi\v2 

book direct ferromagnetic exchange. Note that in transi- 
tion metals / is on the order of 0.5-1 eV. Assuming that 
the tails of the d— wave function decay as exp(— r/r^), 
we see that the ferromagnetic exchange is of the order 
of /exp(— i?/rrf). In t, for weak overlaps, the main role 
is played by the kinetic energy, t ^ — V^|02) ^ 
(01 1^2) /rj ~ exp(-i?/2r<i)/r^. The ratio Jfm/Jafm is 
thus /2exp(-i?/rd)/[exp(-i?/2rd)/r2]2 = /V^. Recall- 
ing that in solids the internuclear distance R is of the 
order of the lattice parameter a, while the Fermi vec- 



l+4>lEst = - P{4>l^V2) 
1 92^ 

'1 ~ <^2)^/4- Thus this term favors ferromagnetism 
2i'^2^ g.j^^ jg ^j^g LSDA countrepart of the text- 
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tor is of the order of the Brillouin zone radius, tt/o, we 
can estimate Jfm / Jafm ~ where 
md ~ 5 — 10 is the effective d— band mass, Ep ~ 3 — 5 eV 
is the d— band width, ~ 1 (in Bohr radii), a ~ 5 — 10. 
Thus Jfm / Jafm ~ 10~^. 
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